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ABSTRACT 

We find that the fraction of classical Broad Absorption Line quasars (BALQSOs) among the FIRST radio 



sources in the Sloan Data Release 3, is 20.5^j \% at the faintest radio powers detected (Li.4ghz ~ 10 Jz erg s '), 
and rapidly drops to < 8% at Li.4ghz ~3 x 10 33 ergs _1 . Similarly, adopting the broader Absorption Index (AI) 



definition of Trump et al. (2006) we find the fraction of radio BALQSOs to be 44^H% reducing to 23.1 
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at high luminosities. While the high fraction at low radio power is consistent with the recent near-IR estimates 
by Dai et al. (2008), the lower fraction at high radio powers is intriguing and confirms previous claims based 
on smaller samples. The trend is independent of the redshift range, the optical and radio flux selection limits, 
or the exact definition of a radio match. We also find that at fixed optical magnitude, the highest bins of 
radio luminosity are preferentially populated by non-BALQSOs, consistent with the overall trend. We do find, 
however, that those quasars identified as AI-BALQSOs but not under the classical definition, do not show a 
significant drop in their fraction as a function of radio power, further supporting independent claims for which 
these sources, characterized by lower equivalent width, may represent an independent class with respect to the 
classical BALQSOs. We find the balnicity index, a measure of the absorption trough in BALQSOs, and the 
mean maximum wind velocity to be roughly constant at all radio powers. We discuss several plausible physical 
models which may explain the observed fast drop in the fraction of the classical BALQSOs with increasing 
radio power, although no one is entirely satisfactory . A strictly evolutionary model for the BALQSO and 
radio emission phases requires a strong fine-tuning to work, while a simple geometric model, although still 
not capable of explaining polar BALQSOs and the paucity of FRII BALQSOs, is statistically successful in 
matching the data if part of the apparent radio luminosity function is due to beamed, non-BALQSOs. 
Subject headings: : black hole physics - galaxies: evolution - galaxies: active - galaxies: jets 



1. INTRODUCTION 

It is a main theme in Active Galactic Nuclei (AGN) stud- 
ies to describe various AGN phenomena within coherent 
schemes. The broad absorption line quasars (BALQSOs) and 
radio quasars (RQs) are both subsamples of quasars, where 
BALQSOs exhibit blue-shifted rest-frame ultra-violet absorp- 
tion troughs (e.g., Weymann et al. 1991) and RQs are char- 
acterized by their radio emission (e.g., Matthews & Sandage 
1963; Schmidt 1963; Kellermann et al. 1989). The fractions 
of BALQSOs and RQs are the subjects of many studies (e.g., 
Weymann et al. 1991; Hewett & Foltz 2003; White et al. 
2007). Recently, Dai et al. (2007, DSS hereafter), adopt- 
ing the classical Balnicity Index (BI) definition (Weymann 
et al. 1991) which selects as BALQSOs those sources with 
absorption troughs wider than 2000kms~ 1 , obtained an in- 
trinsic BALQSO fraction of (20 ±0.2)% using the BALQ- 
SOs detected in the Two Micron All Sky Survey (2MASS) 
and in the Sloan Digital Sky Survey (SDSS) Data Release 3 
(DR3) BALQSO catalog of Trump et al. (2006). DSS showed 
that the BALQSO fraction calculated in the infrared is sig- 
nificantly higher by a factor of ~ 2 than that inferred from 
optical samples alone. Since the infrared bands are less af- 
fected by selection effects due to dust extinction and absorp- 
tion troughs, DSS concluded that the infrared fractions should 
be much closer to the true intrinsic fractions of BALQSOs. 
The intrinsic fraction from DSS is consistent with the early 
estimate by Hewett & Foltz (1993) inferred by directly cor- 
recting the obscuration effects in the optical bands, and subse- 
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quent analysis using similar methods by Knigge et al. (2008). 
In addition, such a a high intrinsic fraction has been recently 
confirmed by Maddox et al. (2008), from a deeper sample 
constructed from the UKIRT Infrared Deep Sky Survey Early 
Data Release, which probes about three magnitudes fainter 
than 2MASS. DSS also showed that an intrinsic BALQSO 
fraction of 43% is found when adopting the broader Absorp- 
tion Index (AI) classification by Trump et al. (2006) which 
identifies BALQSOs as those quasars with absorption troughs 
which have an equivalent width (measured between and 
29000 kms- 1 ) greater than lOOOkms -1 . The validity of the 
latter definition has however been recently questioned (e.g., 
Ganguly et al. 2007; Knigge et al. 2008), and in this paper we 
will therefore discuss results adopting both definitions. 

On the other hand, recent results on the fraction of RQs in 
optical samples are between 10-20%, depending on the def- 
inition (e.g., Sikora et al. 2007 and references therein). In 
several studies, the RQ fraction is found to depend on the 
redshift and optical luminosity (e.g., Jiang et al. 2007). In 
general, constraining the intrinsic fractions of BALQSOs and 
RQs in large, complete quasar samples, is important for un- 
derstanding the nature of these phenomena such as their effi- 
ciency and/or black hole spin distributions (e.g., Murray et al. 
1995; Sikora et al. 2007; Shankar et al. 2008b), and helping 
constrain AGN and galaxy evolutionary models (e.g., King 
2003; Croton et al. 2006; Granato et al. 2006; Hopkins et al. 
2006; Lapi et al. 2006; Shankar et al. 2006). 

The fraction of BALQSOs in RQs is particularly interest- 
ing as it may further constrain which quasar sub-samples are 
caused by geometrical viewing effects and which are caused 
by evolutionary effects, as well as provide additional clues on 
the nature and connections between these two apparently dis- 
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connected quasar phenomena. In a pioneer study, Stocke et al. 
(1992) found no evidence for luminous radio sources among 
68 BALQSOs. In a sample of 255 radio sources in the Large 
Bright QSO Survey, Francis et al. (1993) noted that the frac- 
tion of radio-moderate quasars with BAL features was signifi- 
cantly higher than in the radio-quiet population. Brotherton et 
al. (1998) found five BALQSOs in a complete sample of 111 
QSO candidates also detected in the NRAO VLA Sky Survey. 
Using the FIRST Bright Quasar Survey, Becker et al. (2000) 
found 29 BALQSOs, about 14-18% of their total sample, sug- 
gesting a frequency of BALQSOs significantly greater than 
the typical ~ 10% inferred from optically selected samples at 
that time. Within the small sample, they also found that the 
BALQSO fraction has a complex dependence on radio loud- 
ness. However, the trend they found was driven by the 1 1 low- 
ionization BALQSOs, whereas the 18 high-ionization BALQ- 
SOs were slightly under-represented at high radio loudness, 
but at a statistically insignificant level. With a larger sam- 
ple of 25 high-ionization and 6 low-ionization BALQSOs se- 
lected from an extension of the FIRST Bright Quasar Survey 
to the South Galactic cap and to a fainter optical magnitude 
limit, Becker et al. (2001) confirmed the results by Becker 
et al. (2000) of a steady decrease of the high-ionization 
BALQSO. Using the Large Bright Quasar Survey, Hewett and 
Foltz (2003) argued that with a consistent BALQSO defini- 
tion there is no inconsistency in the BALQSO fractions be- 
tween optically-selected and radio-selected quasar samples. 
They also found that "optically-bright BALQSOs are half as 
likely as non-BALQSOs to be detectable Si.4ghz ^ 1 mJy as 
radio sources". In this paper, we combined the SDSS-DR3 
BALQSO sample (Trump et al. 2006), containing more than 
4,000 BALQSOs, with the FIRST survey data (Becker et al. 
1995) to study the fraction of BALQSOs in the radio detected 
sample. 

2. THE DATA 

The parent sample for the Trump et al. (2006) BALQSO 
catalog is the SDSS DR3 quasar catalog (Schneider et al. 
2005). Data were taken in five broad optical bands (ugriz) 
over about 10,000 deg 2 of the high Galactic latitude sky. The 
majority of quasars were selected for spectroscopic followup 
by SDSS based on their optical photometry. In particular, 
most quasar candidates were selected by their location in the 
low-redshift (z < 3) ugri color cube with its /-magnitude limit 
of 19.1. A second higher redshift griz color cube was also 
used with a fainter /-magnitude limit of 20.2. Following DSS, 
we focus on redshift range ofl.7<z<4.38, where BALQ- 
SOs are identified by C IV absorption in the SDSS spectra. 
In the following we present results using both the Al and BI 
definitions for BALQSOs. 

The DR3 quasar catalog by Schneider et al. (2005) also 
provides the radio flux for those sources which have a coun- 
terpart within 2" in the FIRST catalog (Becker et al. 1995). 
However, we use the integrated flux density observed at 1 .4 
GHz as listed in the FIRST catalog to calculate the specific 
luminosity emitted at 1.4 GHz, Li .4ghz, via the relation 

f. x 1Q-26 

Li.4GHz = 77- u _ a 47rD| ergs-' Hz , (1) 
(i +z) 

where f nt is the integrated flux density in mJy, a is the spec- 
tral index (i.e., f{v) oc v~ a ), and D L is the luminosity dis- 
tance in cm as calculated from the redshift setting 0,„ = 0.3, 
Oa = 0.7 and h = 0.7. We assume a ~ 0.7, which is inter- 
mediate of the typical range for radio jets (Bridle & Perley 



1984). Our results do not qualitatively change if we assume 
either a = 0.5 or a = 1 .0. 

According to Schneider et al. (2005), while only a small 
minority of FIRST-SDSS matches are chance superpositions, 
a significant fraction of the DR3 sources are extended radio 
sources. Therefore the FIRST catalog position for the lat- 
ter may differ by more than 2" from the optical position. To 
perform a comprehensive study of the BALQSO radio frac- 
tion, we built a full FIRST-SDSS cross-correlation catalog, 
containing all the detected radio components within 30" of a 
optical quasar. As our reference, following Schneider et al. 
(2005), we primarily present results for the FIRST-SDSS cata- 
log with radio counterparts defined within 2" , which contains 
a total of 652 radio sources, out of which 78 are BALQSOs, 
about three times the sample used by Becker et al. (2001) 
for BALQSOs. We will discuss results obtained by enlarg- 
ing the radio matches to 5" . In addition, many of the opti- 
cal sources in SDSS are associated to more than one radio 
component in FIRST within 30" , as expected if these sources 
are extended with jets and/or lobes separated from the central 
source. We will therefore discuss results after the inclusion 
of BALQSO sources characterized by multiple radio compo- 
nents, for which we use the sum of the flux densities as a 
proxy for the total radio luminosity of the source. 

Finally, we also stress that FIRST efficiently identifies radio 
matches to optically-selected quasars. By cross-correlating 
SDSS with the large radio NRAO VLA Sky Survey (Condon 
et al. 1998), Jiang et al. (2007) found in fact that only - 6% 
of the matched quasars were not detected by FIRST. Most of 
these sources have offsets more than 5" and are uniformly dis- 
tributed in luminosity thus not altering the main conclusions 
of this paper. 

3. RESULTS 

3.1. The Radio fraction of BALQSOs 

The solid crosses in the left panel of Figure[T]show the frac- 
tion 

, ^Vbal.radio 

/bal = — , (2) 

^VRADIO 

as a function of radio luminosity, where A^radio is the to- 
tal number of RQs detected in SDSS and A^bal.radio is the 
number of RQs which are also classified as BALQSOs. For 
this plot we only used the "complete" SDSS optical sam- 
ple with apparent optical magnitude m, < 19.1. The bins in 
ii.4GHz were created adaptively starting from low radio lumi- 
nosities by requiring a statistically significant sub-sample of 
50 quasars in each bin of luminosity. In each bin we used the 
binomial theorem to estimate /bal along with its ±ler confi- 
dence levels (e.g., Gehrels et al. 1986). 

The first important result of our analysis is the high frac- 
tion of classical (BI) BALQSOs detected at low radio pow- 
ers. The BALQSO fraction of 20^\% at the faintest radio 

powers of Li.4ghz ^ 10 33 ergs~' is in strikingly good agree- 
ment with the 20% classical BALQSO occurrence inferred 
by DSS in the 2MASS-SDSS sample (shown with horizon- 
tal, thin-dotted line in Figures Q] and [3). This result further 
supports and extends the key result discussed in DSS that the 
detected BALQSO fraction should be closer to the intrinsic 
one at longer wavelengths, which are gradually less affected 
by absorption troughs and dust extinction. 

The second striking result of our study is that /bal drops 
rapidly with increasing radio power to 8.OI3 \% at Li.4GHz ^ 
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2 x 10 33 ergs _I . This behavior strongly supports some close 
connection between radio and BALQSO phenomena and, as 
discussed below, poses serious challenges to simple evolu- 
tionary models. 

The open symbols with dashed crosses represent the frac- 
tion /bal obtained including all Al-selected BALQSOs. As 
expected, we get a much larger fraction of BALQSOs within 
our sample of 44^' g%, in very good agreement with the find- 
ings by DSS who found a fraction of 43% AI BALQSOs in 
2MASS (thick-dotted line). Even in this case, /bal presents 
a very similar behavior as a function of radio luminosity and 
that we will discuss in § 13.21 

3.2. Robustness of Results 

Our results are derived by cross-correlating optical quasars, 
RQs, and BALQSOs. Therefore it is natural to expect that 
several issues may affect our conclusions. In this section we 
carefully analyze the most significant sources of bias which 
can enter our computations and find that none of these seems 
to be responsible for the observed /bal — ^i.4GHz behavior. 

We first tested for biases from redshift or luminosity se- 
lection effects. When repeating our analysis using only the 
sources in two redshift bins 1.7 < z < 2.5 and, separately, 
2-5 < z < 4.38, we find very similar trends for the /bal- 
L\ 4ghz relation. 

The completeness of the Schneider et al. (2005) quasar 
sample has been studied in detail by Vanden Berk et al. 
(2005), who estimate a total completeness of 89.3% down 
to m, = 19.1. We do not expect the behavior of /bal to be 
severely influenced by these mild incompleteness effects. As 
extensively discussed (e.g., Cirasuolo et al. 2003; Shankar et 
al. 2008b, and references therein), the radio and optical lu- 
minosities are correlated with a large observed scatter. There- 
fore, incompleteness effects in the optical sample will, on av- 
erage, induce similar effects in the numerator and denomina- 
tor of equation (O, given the extended distribution of optical 
sources at fixed radio luminosity. As an example of the effect 
of different optical luminosity thresholds on the fraction of 
BALQSOs, we show the resulting /bal with no optical mag- 
nitude cut in the right panel of Figure Q] This highly incom- 
plete sample in the optical still shows the exact same drop of 
at least a factor of ~ 4 in /bal- 

On the other hand, our sample may suffer from selection ef- 
fects in FIRST. At 3 mJy the completeness of FIRST is about 
93%, decreasing to 75% at 1.5 mJy, and 55% at 1.1 mJy (see 
Figure 1 in Jiang et al. 2007). To check whether our results 
were affected by this incompleteness, in Figure|2]we compare 
the results shown in the left panel of Figure Q] obtained from 
the full sample (solid crosses), with the sample defined by 
only those radio sources with / nt > 3 mJy (dashed crosses). 
This cut in radio flux density ensures the best completeness 
as possible from FIRST. It is evident that the strong observed 
drop in /bal remains unchanged. 

As anticipated in § [2j the 2" cross-correlation sample en- 
sures a secure identification of the core radio power associ- 
ated to an optical quasar. However, to check for any bias 
in our results due to incomplete radio identifications in the 
SDSS DR3 quasar catalog, we recomputed /bal using the full 
sample built by matching SDSS and FIRST. The left panel 
of Figure [3] shows the fraction of BALQSOs with the radio 
counterpart defined by the nearest source in FIRST within 
5" . The right panel of Figure [3] shows instead the fraction of 
BALQSOs with at least one radio counterpart within 30" of 
the optical source. In both cases we still find /bal to drop 



with increasing radio power. These results are actually ex- 
pected, as a closer inspection revealed that the overall sample 
of radio-optical sources increases by just ~ 4% when includ- 
ing sources within 5", out of which only 3 extra BALQSOs 
are found with respect to our reference sample. The increase 
is even milder (~ 3%) when the 30" sources are included, with 
only 2 extra BALQSOs included in the sample. The paucity 
of extended radio sources in our sample is in agreement with 
Gregg et al. (2006) who found only six quasars from the DR3 
and FIRST Bright Quasar Survey that exhibit Fanaroff-Riley 
Il-type (FRII) morphologies and BAL properties, with a BI 
index strongly decreasing with increasing radio loudness. 

In Figure |4] we show /bal as a function of absolute optical 
magnitude. We do not find significant evidence for any sharp 
trend of /bal as a function of M,-, confirming the results by 
DSS for the 2MASS sample. This implies that different op- 
tical magnitude cuts do not affect /bal on one hand, and, on 
the other, implies that the drop in /bal is mainly a "radio phe- 
nomenon", not linked with the quasar bolometric luminosity. 
Note that we did not apply any correction in Figure |4] for ex- 
tinction as was done in DSS, therefore the resulting fractions 
in the optical are lower than the intrinsic ones (see DSS for 
details). 

The AI sample by Trump et al. (2006) may suffer from a 
significant contamination of false positives, which the close 
inspection by Ganguly et al. (2007) found to be about 30% 
(see also Knigge et al. 2008). To investigate the impact these 
findings may have on our results we have compared the sam- 
ple of sources cataloged as BALQSOs using the AI defini- 
tion, with the sample of BALQSO defined using the classi- 
cal BI index. The distribution of AI-BALQSOs with radio 
counterparts as a function of the AI index (open histogram) is 
compared in the upper left panel of Figure|5]with the distribu- 
tion of the sources which are also identified as BI-BALQSOs 
(filled histogram). We find the latter sample strongly peaked 
at higher AI values, while just the opposite is true for most 
of the AI-BALQSOs. These findings are fully consistent with 
the recent results by Knigge et al. (2008) who claim for a well 
distinct separation for the two BALQSO samples. The upper 
right panel of Figure plots the fraction of AI- but not BI- 
identified BALQSOs as a function of radio power: the drop in 
/bal with increasing radio power is much less evident for this 
subclass of sources compared to BI-BALQSOs. These results 
suggest that AI-BALQSOs are a different class of quasars 
with respect to the BI-BALQSOs. In the lower right panel 
of Figure [5] the dashed, dotted and solid lines show the frac- 
tions of AI-BALQSOs, AI- but not BI-BALQSOs, and BI- 
BALQSOs, respectively. This plot confirms that the apparent 
drop in the AI-BALQSOs with radio power is driven by the 
behavior in BI-BALQSOs. The lower right panel of Figure[5] 
shows the cumulative distributions for the different subclasses 
of BAL and non-BALQSOs, as labeled. We find that the prob- 
ability that AI but not BI selected BALQSOs are drawn from 
the same distribution as non-BALQSOs is P = 0.26, while the 
probability for being drawn from any other distribution drops 
by several orders of magnitude. 

In Figure|6]we show the distributions of all non-BALQSOs 
(dashed lines), BI-BALQSOs (solid lines), and AI- but not BI- 
BALQSOs (dotted lines) as a function of radio luminosity for 
different bins of absolute optical magnitude, as labeled. These 
bins were chosen to provide an equal number of BALQSOs 
in each panel. It is apparent that compared to BI-BALQSOs, 
both the AI- but not BI-BALQSOs and the non-BALQSO dis- 
tributions have a greater tendency to populate significantly 
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higher, potentially beamed, radio luminosity bins. This pro- 
vides further evidence that the trends in /bal are not a mere 
consequence of random effects in the radio power distribu- 
tions of optical quasars. The dashed, dotted, and solid lines 
in Figure [7] show, respectively, the histograms of the distribu- 
tions of BI-BALQSOs, AI- but not BI- BALQSOs, and non- 
BALQSOs as a function of color defined as the difference 
between absolute /-band and AB radio magnitudes listed in 
Schneider et al. (2005). This plot shows that on average AI- 
but not BI- BALQSOs, and non-BALQSOs have very similar 
optical-to-radio distributions while BI-BALQSOs are clearly 
offset populating lower radio luminosity bins at fixed optical 
magnitude. 

Overall our results show that the classical BALQSO frac- 
tion in the radio is as high as inferred from the near-IR and 
rapidly drops from ~ 20% to < 4% with increasing radio 
power. This clear trend confirms previous claims by Becker 
et al. (2001), taking advantage of a sample of ~ 3 higher. 
We conclude that the behavior of /bal as a function of ra- 
dio power seems to be an intrinsic, real feature of the radio, 
BI-BALQSO population, not linked to any significant selec- 
tion effect. In the following, if not alternatively specified, we 
will always refer to BALQSOs as those defined through the 
classical BI definition. 

4. EXPLAINING THE OBSERVED TREND 

Models aimed at explaining the BALQSO and/or radio frac- 
tions, can be broadly divided into two classes: one in which 
each quasar phase is distinct and characterized by the inter- 
nal physical evolution undergone by the system, and the other 
in which all these phenomena co-exist at all times but differ- 
ent source orientations modify the observed balance of one 
process with respect to another. 

4. 1 . The Evolutionary Model 

The results discussed in § [3] seem to pose a serious chal- 
lenge to "evolutionary" models for BALQSOs. The simple 
sketch shown in the left panel of Figure [8] represents the ba- 
sic picture for the evolution of a black hole within its host 
galaxy (e.g., Kawakatu et al. 2003). The mass of the central 
black hole, shown as a solid line in the left panel of Figure [8] 
exponentially grows through gas accretion from the surround- 
ing medium. Initially the whole galaxy, or at least its central 
region, is buried under optically-thick layers of dust which 
allow the AGN to be detectable only in those bands that suf- 
fer less from obscuration, such as in X-rays (e.g., Ueda et al. 
2003). When the black hole becomes sufficiently luminous, 
it injects enough energy and momentum (e.g., Fabian 1999) 
into the surrounding medium to disperse the obscuring mate- 
rial. During this blowout phase, the system will be detected as 
a BALQSO featuring the signatures of strong winds in their 
spectra, while at later stages it will be observed as an optical, 
dust-free quasar. 

It is still unclear what triggers the radio activity of an AGN 
(e.g., Sikora et al. 2007), so it is difficult to predict when 
an optical quasar should become radio-loud (e.g., Blundell & 
Kuncic 2007; Shankar et al. 2008b). Nevertheless, in pure 
evolutionary models the radio phase is usually viewed as a 
brief period within the optical phase as the incidence of RQs 
in optical samples is rather small, (~ 15 — 20%, e.g., Jiang et 
al. 2007; see also § [TJ. In this model, our results on /bal 
indicate that the radio duty cycle must overlap the boundary 
between the BALQSO and optical phases, because a signifi- 



cant fraction of optical quasars also emit in radio but are not 
BALQSOs, and vice versa. 

In a pure evolutionary model our result that /bal ~ 20% 
constrains the BALQSO and non-BALQSO phases to last for 
~ 20% and 80%, respectively, within the brief radio phase, 
which is about 15-20% of the total optical quasar lifetime (see 
left panel of Figure [8). The consistency between the optical 
and radio BALQSO fractions imposes a coincidence problem, 
which could be otherwise naturally interpreted in orientation 
models (§ I4.21 i. Additionally, the drop in /bal at high radio 
powers is difficult to reproduce simply in pure evolutionary 
models, and thus such models would need to be fine-tuned to 
explain the drop. One simple explanation for such a trend may 
be an increasing feedback efficiency with quasar luminosity 
(e.g., Kelly et al. 2008). The more luminous sources would 
then be faster in removing the surrounding obscuring medium, 
shortening the BALQSO phase and consequently inducing a 
decrease in /bal- However, this model is at variance with the 
results in Figure|4]and by DSS which show that the BALQSO 
fraction is constant with optical luminosity. 

Gregg et al. (2006) found an absence of BALQSO FRII 
sources and a drop in the BI index with increasing radio loud- 
ness of the sources. They interpreted their results within an 
evolutionary scheme where a BALQSO occurs in the early 
stage of a quasar during its emergence from a thick shroud of 
surrounding dust (see left panel of Figure [8]). In this model 
FRIIs are only seen after most of the obscuring shroud has 
been removed. Following Gregg et al. (2006) we further 
probe the constraints for evolutionary models by address- 
ing the coexistence of FRII/BAL sources within our sample. 
In Figure [9] the open triangles show the fraction of FRII- 
selected sources 2 within the whole radio sample, while the 
filled circles refer to the fraction of FRII sources within the 
Bl-selected radio sample. Both samples have been divided 
into two subgroups with radio luminosities below and above 
^i.4GHz = 2 x 10 33 ergs -1 Hz -1 . Note that the contamination 
from random interloping radio sources in the identification 
of FRII sources is small. The average number of interlopers 
within 30" of any position is only 0.02 (Becker et al. 1995), 
i.e., a cumulative number of just about 3 sources are expected 
to be misidentified as FRIIs. It can be seen that in both BI 
and not Bl-selected samples, the fraction of FRII increases 
with radio power, mirroring the fact that FRIIs are intrinsi- 
cally more luminous radio sources (e.g., De Zotti et al. 2005). 
We do find some evidence for FRIIs to be less associated to 
BALQSOs within the Bl-selected sources than in the overall 
sample at fixed radio power, although the difference is only 
marginally significant. Gregg et al. (2006) claim a drop by a 
about an order of magnitude in the fraction of FRIIs in BAL- 
radio samples, while we find a mean drop of ~ 2.4. Also, we 
do not find strong evidence for the Gregg et al. (2006) anti- 
correlation between BI index and radio loudness. A Pearson's 
r-test yields a correlation parameter of r = —0.36, with a sig- 
nificance of P — 0.55; however, the small number of FRII 
BALQSOs limits this analysis. The BI index of FRIIs in our 
BAL-radio sample does not show any clear trend with radio 
power, as shown with filled circles in the right panel of Fig- 
ure [9] For comparison, in the same Figure we also show with 
open symbols the mean Bis for the low and high radio power 
subsamples of non-FRIIs, which show a rather flat behavior 
with increasing radio power, in agreement with the overall 

2 Here we simply call "FRIIs" all those sources in the sample with multiple 
radio components within 30" . 
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sample (see § !4.3l and Figure [TTTi. 

It also may well be true that the BALQSOs we are referring 
to here are only a partial representation of the "outflowing" 
AGNs. Many more systems may be detected during initial 
and final phases of the blowout and therefore may be char- 
acterized by lower velocities. Ganguly & Brotherton (2008) 
have shown that including Narrow (< 800km s _1 ) and Asso- 
ciated absorption systems, the latter being characterized by 
narrow absorption-lines near the redshift of the quasar, the 
fraction of outflowing AGNs rises to about 60%, independent 
of optical luminosity. The /bal radio dependence for the cu- 
mulative sample may then be different then what is inferred 
here. Several works have shown that a significant fraction of 
RQs show associated absorption and high radio power sources 
may preferentially exhibit narrower absorption features (e.g., 
Weymann et al. 1979; Anderson et al. 1987; Ganguly et al. 
2001 ; Vestergaard 2003). Such a hypothesis needs to be tested 
with samples with identified Associated Absorption systems, 
such as those by Ganguly et al. (2007). Nevertheless, the drop 
in /bal for BI-BALQSOs will still need to be explained. 

4.2. The Orientation-Beaming Model 

In this section we show that a simple orientation effect can 
fully explain the puzzling trend in the /bal-Li.4GHz distribu- 
tion, although even this model is not entirely satisfactory, as 
discussed at the end of the section. This model assumes that 
some fraction of RQs are relativistically beamed towards the 
observer and are boosted to higher radio luminosities. Since 
the radio luminosity function is steep, beamed radio sources 
are a disproportionate fraction of the bright sources. Accord- 
ing to popular disk wind models (e.g., Murray et al. 1995, 
Proga et al. 2000), radio BALQSOs are instead on average 
non-beamed sources, being preferentially viewed close to the 
plane of the accretion disk surrounding the central black hole. 
It is then reasonable to expect that if BALQSOs are a fixed 
fraction of the intrinsic radio luminosity function, then the 
observed occurrence of BALQSOs must decrease at high ra- 
dio powers due to the apparent increase of A^radio in Eq. (f2]i 
towards high luminosities due to beaming. 

We quantitatively explore this basic idea by building a sim- 
ple "unification" model (e.g., Urry & Padovani 1995). The 
intrinsic radio quasar luminosity function was assumed to be 
a power-law oc with arbitrary normalization and 

slope (3 = —3.0, which is consistent with the bright-end slopes 
measured for the optical and radio quasar luminosity func- 
tions (e.g., De Zotti et al. 2005, Richards et al. 2006, Jiang 
et al. 2006) 3 . We then assume that the population of radio 
quasars is composed of three sub-populations, as sketched 
in the right panel of Figure [8] The first one is represented 
by radio BALQSOs, which are preferentially viewed at lines 
of sight close to the disk or the torus, at a maximum angle 
of 9 < 11.5°, corresponding to a fraction of sin(f?) = 20%, 
as inferred in §[3] and by DSS. Following Urry & Padovani 
(1995) we then assumed that the fraction of beamed sources 
is viewed at an angle of 38° from the jet axis, corresponding 
to a fraction of ~ 21%. The last sub-group is composed of 
radio non-beamed sources viewed at intermediate angles ac- 
counting for the remaining ~ 59% of the total. AI-BALQSOs 
in this model must be randomly distributed at all angles to 
ensure their fraction to be flat with radio power. The relative 

3 The exact value of this slope does not affect our results as the uncer- 
tainties in this parameter are degenerate with the uncertainties in other model 
parameters 



contributions of the different type of sources to the total radio 
luminosity function are shown in the left panel of Figure [TOl 
plotted as L 3 $(L) which flattens out the luminosity function 
and enhances the difference among the lines. 

The luminosity function of beamed sources was computed 
following the method outlined in Urry & Padovani (1991). 
The observed luminosity function at a given observed lumi- 
nosity L bs was obtained by convolving the intrinsic luminos- 
ity function with the probability P(L b s |Lint)- Following Urry 

6 Padovani (1991) we allowed beamed sources to also have 
a diffuse component from the lobes and the core, in addition 
to the beamed emission. Therefore the observed luminosity 
of this sub-group was parameterized as L b s = (1 + f8 p )L mt , 
where L mt is the intrinsic luminosity, / is the fraction of lu- 
minosity in the jets, 5 — [7(1 — /Jcosf?)] -1 is the Doppler 
factor corresponding to a given angle 8 and Lorentz factor 

7 = (1 — /3 2 )~'/ 2 , and the exponent p is a parameter depend- 
ing on the spectrum and reacceleration of the radiating parti- 
cles. Given the relatively high number of parameters in this 
model, we adopt p = 3, 7 = 11,/ = 0.005 and the maximum 
beaming angle to 6* max = 38°. These are the same values used 
by Urry & Padovani (1995) and, more recently, by Padovani 
et al. (2007) to reproduce within unification models the lumi- 
nosity distributions of samples dominated by beamed sources. 

The apparent luminosity function of beamed radio sources, 
"crosses" the intrinsic luminosity function at a luminosity 
which depends on the minimum luminosity in which we 
allow beamed sources to appear in the luminosity func- 
tion (see Figure ITOb . We parameterized this luminosity as 
logLi.4GHz,MiN =32- AlogL L4 GHz, where logL L4G Hz ~ 32 
is the minimum luminosity in the sample. The solid line in 
the right panel of Figure [To] shows our best-fit model for the 
expected fraction of BALQSOs as a function of radio lumi- 
nosity. Optimizing the model to the L\ 4GHz-dependence of 

/bal, we find AlogL 1 . 4G Hz = 0.09^ 2 , with X L ~ 12.3, 
for 12 degrees of freedom. This model implies that beamed 
sources are almost absent in the luminosity function below 
Li.4GHz~ 10 32 ergs~' and then start to become important only 
above this luminosity threshold. To better compare with the 
data, we also show with open circles the expected fraction for 
this model integrated over the same radio luminosity bins of 
the data. 

In this model the dependence of the BALQSO fraction on 
radio luminosity is entirely due to geometry and it is inde- 
pendent of the duty cycle of RQs within the overall quasar 
population. Nevertheless, the origin of the radio emission can 
still be an evolutionary phase or related to specific properties 
of quasars such as the magnetic field, the black hole mass, 
spin or Eddington ratio. 

A pure geometrical model, although plausible on statistical 
grounds, is not entirely satisfactory. Polar BI-BALQSO out- 
flows have indeed been observed in several cases (e.g., Zhou 
et al. 2006; Ghosh & Punsly 2007; Wang et al. 2008), and the 
paucity of FRIIs within BALQSOs can hardly be accounted 
for in simple orientation models (see § 14.11 ). Nevertheless, 
we expect that some non-orientation processes may effect the 
physics of radio and BALQSOs and induce different evolu- 
tionary trends in different type of sources. For example, it 
has been observed that AGNs have broad Eddington ratio dis- 
tributions at fixed black hole mass (e.g., Shen et al. 2008). 
The more luminous FRII sources may be driven by higher Ed- 
dington ratios than FRI sources, at fixed BH mass, and con- 
sequently may have higher kinetic powers and stronger radia- 
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tion pressures. These may be more efficient in removing BAL 
clouds along the line of sight and also be more efficient in 
driving collimated jets, thus explaining both the reduced BAL 
fraction and the large scale hot-spots observed at kpc scales 
(as also partly discussed by Gregg et al. 2006). 

4.3. The Energy Exchange Model 

Another possible explanation for the drop in /bal, could 
rely on some sort of "energy exchange" from the wind to the 
jet in high luminosity radio sources. Within the framework of 
the geometrical model, the energy gained by the system from 
the external gas accretion must be redistributed between the 
jet and the wind by some physical mechanism, such as the 
magnetic field. We could then expect that those quasars with 
high radio powers preferentially inject energy into the jet due 
to a stronger intrinsic magnetic field and/or a rapidly spinning 
central black hole, thus limiting the energy channeled into the 
wind and inevitably causing a drop in /bal- 

We tested this idea by considering in the left panel of Fig- 
ure[TT]the BALQSO BI index (upper panel), AI index (middle 
panel) and maximum velocity V max (lower panel) measured 
from the absorption troughs as proxy of the wind power, and 
studied their behavior as a function of radio luminosity. In 
each bin, the quantities were computed from the biweight- 
mean (Hoaglin et al. 1983) and the errors on this mean were 
estimated from the biweight-er reduced by V^V, where N is the 
number of BALQSOs in the bin. Becker et al. (2000, 2001) 
found the mean BI decreases with radio luminosity for the 
high-ionization BALQSOs (HIBALQSOs), an effect which 
could be simply interpreted as a decrease in the amount of 
obscuring mass for the more luminous sources. They also 
found tentative evidence for HIBALQSOs to have their maxi- 
mum velocity decreasing with radio power. In our larger sam- 
ple, which should be dominated by HIBALQSOs, we instead 
find the AI, BI and maximum wind velocity show no clear 
evidence for a decrease at high radio powers although with 
a large scatter, in conflict with the energy exchange model. 
We also find marginal evidence for a drop in the AI index but 
not in the V max of the AI but not Bl-defined BALQSOs (right 
panel of Figure [TT). 

The energy exchange model may also conflict with evi- 
dence that AGNs have similar radiative and kinetic efficien- 
cies (e.g., Shankar et al. 2008a,b). Moreover, if the efficiency 
is constant, then V max should increase with luminosity (i.e., 
with the mean accretion rate), rendering the result on its flat- 
ness rather intriguing (e.g., Ganguly et al. 2007). 

5. DISCUSSIONS 

Synchrotron self-absorption may impose an uncertainty in 
estimating the radio flux density. The spectral analysis by 
Becker et al. (2000) found that about one third of the sources 
could show synchrotron self-absorption at lower frequen- 
cies, typical of compact radio sources. If so, up to 30% 
of the radio luminosities of BALQSOs may be underesti- 
mated and in turn significantly affect our results. To quan- 
tify the implications of self-absorption, one could assume 
the intrinsic fraction of BALQSOs to be /bal = k, with k 
constant with increasing radio luminosity. If self-absorption 
mainly affects sources above an intrinsic radio luminosity 
of Lbreak = £i.4GHz ~ 2 x 10 33 ergs _1 , as suggested by the 
observed behavior in /bal, then the observed fraction of 
BALQSOs, /sALobs, will be "distorted" by the fact that a 
fraction of ~ 30% of luminous radio sources will be trans- 



ferred from high to low luminosities, decreasing the num- 
ber of BALQSOs at high radio powers and proportionally 
increasing it at low radio powers. For example, by setting 
k = 0.15, the fraction of BALQSOs with L\ 4ghz ^ ^break will 
be /bal obs = [k x (1 - 03)]/[k x (1 - 0.3) + (1 - k)] ~ 0.1, 
while the fraction of sources with L\ 4ghz < ^break will be 
/bal obs = [*x (l + 0.3)]/[Jtx (1 + 0.3) + (1-*)] ~ 0.19, 
which is only in marginal agreement with observations. This 
simple explanation therefore cannot fully account for the 
steep drop in /bal obs- To be relevant in our study, synchrotron 
self-absorption should also be able to decrease the intrinsic 
radio power by factors of 10 to 1000, however typical reduc- 
tions of only a factor of ~ 10 are expected (e.g., Snellen & 
Schilizzi 1999). 

6. CONCLUSIONS 

We find that the fraction of classical BALQSOs among 
the FIRST radio sources in the SDSS DR3 catalogue, is 
20.5^5 9% at the faintest radio powers detected (Li.4GHz ~ 
10 32 ergs~'), and rapidly drops to < 8% at Li.4ghz ~ 3 x 
10 33 ergs _1 . Similarly, adopting the broader AI definition of 
Trump et al. (2006) we find the fraction of radio BALQSOs 
to be 44^7 g% reducing to 23.llg 3 at high luminosities. 

The detected fraction at low radio luminosities is consistent 
with the recent estimates inferred in infrared bands by DSS, 
supporting the fact that longer, less absorbed wavelengths, are 
more suitable for probing the intrinsic fraction of BALQSOs. 
We therefore agree with Hewett and Foltz (2003) who claimed 
a similar occurrence of BALQSOs in optical and radio sam- 
ples. The variation of the BALQSO fraction with radio power 
is also similar to that found by Becker et al. (2000, 2001) and 
Hewett and Foltz (2003), although our results are supported 
by a much larger sample. The decrease in the number of ra- 
dio BALQSOs is real, in the sense that it is not dependent on 
redshift or luminosity cuts, biased by optical or radio selec- 
tion effects, or dependent on the specific type of radio source 
considered (core-dominated or extended with multiple com- 
ponents). We also find that at fixed optical magnitude, the 
highest bins of radio luminosity are preferentially populated 
by non-BALQSOs, consistent with the overall trend. 

However, we find that those quasars identified as AI- 
BALQSOs but not under the classical BI definition, do not 
show the same significant drop in the fraction as a function 
of radio power and share similar cumulative radio distribu- 
tions as non-BALQSOs. This further supports independent 
claims for which these sources characterized by lower equiv- 
alent width, may represent an independent class with respect 
to the classical BALQSOs. We find the BI, AI and mean max- 
imum wind velocity to be roughly constant at all radio powers. 

We discuss several plausible physical models which may 
explain the observed fast drop in the fraction of the clas- 
sical BALQSOs with increasing radio power, although no 
one model is entirely satisfactory. A strictly evolutionary 
model for the BALQSO and radio emission phases requires 
a strong fine-tuning to work, while a simple geometric model 
where the apparent radio luminosity function is partly due to 
beamed, non-BALQSOs, although more promising in match- 
ing the data, cannot explain strong polar BALQSOs and the 
paucity of FRII sources within BALQSOs. Of course, the 
drop in /bal could also be viewed as an increase of the 
BALQSO fraction towards low-intermediate radio sources 
which may be more likely to show BALQSO features. This 
model however does not explain why the fraction of BALQ- 
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SOs at low radio powers is identical to the one obtained in the 
optical-NIR. 
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proved the paper. We also thank Chris Kochanek for helpful 
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sion at the AGN lunch in the OSU astronomy department. 
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FIG. 1 . — Left panel: fraction of BALQSOs within the radio sample of SDSS quasars with optical luminosity brighter than m, = 19.1 as a function of radio 
luminosity and with radio counterpart within 2" of the optical source; the filled symbols with solid lines show the sample of BALQSOs defined under the BI 
definition, while the open symbols with dashed lines show the sample of BALQSOs sample compiled following the Al definition; the lower dotted and upper 
dotted lines refer to the "intrinsic fractions" of 20% and 43% of BALQSOs found in 2MASS by Dai et al. (2007), respectively. Right panel: fraction of BALQSOs 
within the full SDSS sample of RQs; here all lines and symbols are as in the left panel; the statistical significance of the drop in the radio-BALQSOs fraction 
towards high radio powers does not depend on the optical magnitude limit. 
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FIG. 2. — Dependence of /bal on radio luminosity for samples defined to have different flux density limits f mt in FIRST and corresponding different 
completeness levels. The drop in /bal is always present. 
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FIG. 3. — Left panel: fraction as a function of radio power of radio-BALQSOs with an optical magnitude brighter than m, = 19.1 with a single radio spot 
within 5" of the optical source; if multiple radio components are contained within this aperture, then only the closest source is considered. Right panel: fraction 
of radio-BALQSOs for sources which lie within 30" of the optical source; if multiple radio components are present within this aperture, then they are treated as 
a single source with radio power given by the sum of the integrated flux densities of each component. In both panels all symbols and lines have same meaning as 
in Figure[T] 
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FIG. 4. — Fraction of radio-BALQSOs as a function of AB optical magnitude, /bal is about constant if not increasing at higher optical luminosities, which does 
not supports any direct obvious link in the drop of /bal with increasing radio power. Note that the true fractions in the optical should be higher once corrected 
for absorption and extinction (Dai et al. 2008). 
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FIG. 5. — Upper left panel: histogram of all BALQSOs defined with the Al definition; the striped histogram refers to the fraction of the sources which are 
identified as BALQSOs also under the BI definition. Upper right panel: fraction of sources defined as BALQSOs with the Al definition but not with the BI one as 
a function of radio power; the drop in /bal with increasing radio power is much less evident for this subclass of sources. Lower left panel: the dashed, dotted and 
solid lines show the fractions of radio-BALQSOs defined with the Al definition only, with Al but not BI, and with the BI definition, respectively, as a function of 
radio power. Lower right panel: cumulative distributions for the different subclasses of BAL and non-BALQSOs, as labeled, as a function of radio power. The 
probability that Al but not BI selected BALQSOs are drawn from the same distribution as non-BALQSOs is P = 0.26. 
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FIG. 6. — The dashed, dotted, and solid lines show, respectively, the histograms of the distributions of BI-BALQSOs, Al but not BI BALQSOs, and non- 
BALQSOs as a function of radio luminosity in different bins of optical luminosity, as labeled. Unlike the BI BALQSOs, the non-BI BALQSOs have a greater 
tendency to form a prolonged tail of very luminous, potentially beamed, radio sources. 
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FIG. 7. — The dashed, dotted, and solid lines show, respectively, the cumulative histograms of the distributions of BI-BALQSOs, AI but not Bl BALQSOs, 
and non-BALQSOs as a function of color. This plot shows that on average at fixed optical luminosity the non-BI BALQSOs have a tendency to be more radio 
luminous than Bl BALQSOs. 




FIG. 8. — Left panel: schematic diagram of a pure "evolutionary" model. The mass of the central black hole starts growing in a heavily obscured phase. When 
it becomes sufficiently luminous to blow away the surrounding medium, it can then be detected as a BALQSO and then finally as an optical, dust-free quasar. 
Right panel: schematic representation of a simple "orientation" model for RQs and BALQSOs. The radio phase in both schemes is always considered as a brief 
period within the optical phase. 
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FIG. 9. — Left panel: fraction of FRII-selected sources among the Bl-selected BALQSOs (filled circles) and non-BALQSOs (open triangles); both samples have 
been divided into two subgroups with radio luminosities below and above L\ 4ghz = 2 x 10 33 ergs -1 Hz - 1 ; we display the median luminosity of each subsample. 
Right panel: absorption index for non-FRII radio sources (open circles) and radio-BI selected FRII sources (filled circles) a s a f unction of radio power. For the 
more numerous non-FRIIs, we display the mean Bl values and associated uncertainties following Hoaglin et al. (1983; see § 14.31 . 
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FIG. 10. — Left panel: radio quasar luminosity functions plotted in units of L 3( 1>(L) to enhance the differences between the lines; the long-dashed line is 
the luminosity function of RQs detected as BI BALQSOs, normalized to be 20% the one of the intrinsic total radio quasar luminosity function, which we have 
assumed to have a slope of —3.0; the dotted line is the fraction of ~ 59% of radio non-beamed, non-BI BALQSOs, while the dot-dashed line is the fraction 
of ~ 21% of beamed sources (see text for details); the solid line is the total apparent luminosity function. Right panel: the solid line is the expected fraction 
of BALQSOs as a function of radio luminosity derived as the ratio of the dashed and solid lines in the left panel; the open circles show the expected fraction 
averaged over the same radio luminosity bin as the data; all other symbols and lines are as in the left panel of Figure[T] 
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FIG. 1 1 . — Left panel: the crosses in the upper panel show the mean BALQSO Balnicity Index per bin of radio luminosity, while the full sample is shown with 
points; the middle panel shows the same for the Absorption Index, while the lower panel shows the maximum velocity measured from the absorption troughs as 
a function of radio luminosity. Right panel: Absorption Index (upper panel) and maximum velocity (lower panel) for AI but not Bl-selected BALSQSOs. We do 
not find any significant trend with radio luminosity for all plotted quantities. 



